A super-Ohmic energy absorption in driven quantum chaotic systems 
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We consider energy absorption by driven chaotic systems of the symplectic symmetry class. Ac- 
cording to our analytical perturbative calculation, at the initial stage of evolution the energy growth 
with time can be faster than linear. This appears to be an analog of weak anti-localization in dis- 
ordered systems with spin-orbit interaction. Our analytical result is also confirmed by numerical 
calculations for the symplectic quantum kicked rotor. 
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Introduction. — The problem of energy absorption in 
a system driven by an external time-dependent field is 
fundamental and important in many areas of modern 
physics. For a metallic sample of the volume V in an 
external electric field E{t) = Eocosujt the textbook so- 
lution to this problem is given by the expression for the 
Joule heating in the Ohmic regime: Wq = V ao Eq/2, is 
the constant energy absorption rate determined by the 
Drude conductivity ctq- 

This classical picture is based on the linear response 
theory for systems with essentially continuous spectrum 
of electron states. For quantum systems with few de- 
grees of freedom or for mesoscopic systems with many 
degrees of freedom but still appreciable level separation 
the Ohmic regime may break down leading to a time- 
dependent absorption rate W{t). 

This point can be illustrated by an example of the 
quantum kicked rotor (QKR) with the Hamiltonian: 
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KV{e) ^ 5{t-nT), V{e)^cose, (1) 



where (. = —id/dO is the angular momentum, / is the 
moment of inertia, and K is a, constant controlling the 
strength of perturbation. For generic sufficiently large K 
the classical dynamics described by the Hamiltonian (^ 
is completely chaotic. The period-averaged energy ab- 
sorption rate Wq = K"^ /{^TI) in this case is indepen- 
dent of time, analogously to the Ohmic absorption. Yet 
at sufficiently long times t ^ t^, ^ K^P/T the Ohmic 
regime breaks down because of the accumulation of quan- 
tum corrections and the absorption rate decreases to zero. 
This effect is known as dynamic localization (DL) in the 
energy space Jj, and is analogous to Anderson localiza- 
tion for disordered systems |3| ■ Such behavior is not spe- 
cific to QKR, it occurs in other chaotic systems [3, ^ . 

However, if T/(47r/) takes a rational value, the separa- 
tion between certain energy levels of the rotor becomes an 
integer multiple of the frequency 27r/T, and the absorp- 
tion rate is linear in time: W{t) ex t ij. The same takes 
place for a harmonic oscillator coupled to the external 
harmonic field via the coordinate, when the frequency 



of the field is exactly at resonance with the oscillator 
frequency |5|. Fermi accelerator is another example of 
a system where W{t) can grow with time [6|. Such an 
anomalous {growing with time) super- Ohmic behavior is 
typical of resonances. In contrast to the localization, it 
is analogous to the ballistic transport through resonant 
levels in a tight-binding model of Id crystals. One can 
trace it back to the classical integrability of the system 
with a time-dependent perturbation. 

In this letter we consider a class of chaotic systems 
without resonances which show the super-Ohmic energy 
absorption. Namely, we focus on the quantum, correc- 
tions SW{t) to the energy absorption rate in the time- 
dependent random matrix theory (RMT) of the symplec- 
tic symmetry class, described by the Hamiltonian: 



H(t)^Ho + V4>{t), 



(2) 



which possesses the time-reversal symmetry, but not the 
spin-rotation symmetry. Here Hq and V are random ma- 
trices J3l whose symmetry will be specified below, and 
(f>{t) is a given function of time. This model describes 
e. g. the dynamics of electrons in driven quantum dots 
in the presence of a spin-orbit interaction. 

The corresponding problem for the orthogonal and uni- 
tary symmetry classes has been recently considered |8j 
and ana/ytica/ expressions for SW{t) have been obtained. 
For the harmonic perturbation (pit) = cos tot switched on 
at i = 0, the absorption rate W{t) appears to be related 
to the frequency-dependent diffusion coefficient D(uj) in a 
quasi- Id disordered wire of the corresponding symmetry 
class, with the quantum corrections included: 



W{t) 

Wn 



+°" du^ 
.^ 27r 



0) Do ' 



(3) 



where Dq is the classical diffusion coefficient. This re- 
lationship does not contain any specific feature of the 
model and is also valid for the QKR in the region of pa- 
rameters where it can be mapped onto the quasi-ld non- 
linear cr-model ^]. If Eq. Q is valid in the symplectic 
case as well, the energy absorption rate W{t) should grow 
with time beyond the Ohmic limit Wo, as D(lu) is known 



to have positive quantum corrections in the presence of 
a spin-orbit interaction [lOJ. Our calculations presented 
below show that this is indeed the case. 

Choice of the model. — We adopt the following single- 
electron Hamiltonians, which turns out to be the most 
convenient technically for perturbative calculations: 



-J 

Hit) ^^ + U{r) + UUr) + V{f) ^(t), 



(4) 



where p = — iV, and {/(r) and V{r) are independent 
Gaussian random fields: {U{r)U{r')) = au5[f — r'), 
{V{r)V{r')) = avS{f — r'). The spin-orbit interaction 
is also taken to be random [ifl uM '■ 



Usoir) = S V[/so(r) xp 



(5) 



where (Usoir) Uso{r')) = asoSif-r), and a = {a'^,cry,a^) 
are Pauli matrices. 

The advantage of the model Q) is that the spin-orbit 
coupling and the driving perturbation are represented 
by random locally correlated fields. This makes it pos- 
sible, after a proper re-formulation |l2l Il3l | within the 
Keldysh formalism [lj| , to apply basic rules of the impu- 
rity diagrammatic technique [15| and its extension used 
in the theory of weak Anderson localization and meso- 
scopic phenomena [la| to consider essentially nonlinear 
in the driving perturbation, non-equilibrium problems. 

In the absence of the time-dependent term Eq. Q is 
a basic model for describing disordered metals with a 
random spin-orbit interaction. The kinetic energy term 
determines the bulk density of states v (per unit vol- 
ume, per spin projection). Then ajj — l/(27rj/To) and 
[pxp']'^aso = i/i'2'^i^Tso), where tq is the momentum 
relaxation time, Tso 3> tq is the spin relaxation time, 
and [p X j?]"^ denotes the momentum product averaged 
over the Fermi surface. It has been shown [l2| that 
for a finite sample in the long-time, low-energy limit 
e <C 1/tso <C Eth [Erh being the Thouless energy) this 
model reduces to the zero-dimensional nonlinear cr-model 
which, in turn, is equivalent to the RMT of the symplec- 
tic symmetry class. This corresponds to Hq in Eq. (0 
being a random matrix from the Gaussian symplectic 
ensemble (GSE) with the mean energy level separation 
5^ l/(2i/V). 

The equivalence of the model Q to the time-depen- 
dent RMT can be also demonstrated, but it requires ad- 
ditional conditions: 



— = 4rro < -^ 



<1. 



(6) 



Here we have introduced a parameter F = ixvay jl. We 
will always be interested in the limit 5 <g; F <C w, in which 
case 1/F is the time required to absorb one photon, as 
given by the Fermi Golden Rule. The first condition © 
allows to neglect the multi-photon absorption processes 





FIG. 1: DiagramiTiatic representations for (a) diffuson, (b) 
cooperon; the Greek indices label the spin projections f , J,. 



which can easily violate the condition of small energy 
transfer Ae <C Etk even when the condition lo <^ E^h is 
fulfilled. Under the conditions lO the last term of Eq. Q 
corresponds to V<\)it) of Eq. Q with V being a random 
matrix from the Gaussian orthogonal ensemble (GOE), 
whose matrix element mean square is given by FJ/tt. 

We note that if instead of the time-dependent pertur- 
bation with random V(f) one considers a more physi- 
cal form of the perturbation with some fixed V{f^ corre- 
sponding to a uniform electric field or a modulation of the 
quantum dot confinement potential, one can show |l2lll3| 
that under the conditions (|SJ), where only the basic sym- 
metries (time-reversal, spin-rotation) matter, it will be 
equivalent to the same RMT with a properly redefined F. 
However, calculations would be quite cumbersome in this 
case because of the necessity to take into account bound- 
ary conditions. The coupling by a random field V{r) with 
local correlations helps to avoid this technical problem. 

We also note that instead of Eq. (0J one can con- 
sider a model where spin-orbit interaction is in the time- 
dependent term: H{t) = j) /2m + U{r) + Usoir) Ht). 
It corresponds to a different kind of the time-dependent 
RMT, where in Eq. 10) Ho is taken from GOE and V from 
GSE. Still, the result turns out to be the same as for the 
model (0}, being remarkably robust and independent of 
whether the spin-rotational invariance is broken in the 
time-independent or in the time-dependent term [see also 
the discussion after Eq. I|12|) ]. 

Weak dynamic antilocalization. — The model (0J al- 
lows for a perturbative treatment which is very similar 
to the theory of weak Anderson localization Igj. The 
building blocks of this theory are ladder diagrams: the 
diffuson and the cooperon (Fig. QJ. The latter can also 
be represented as the maximally crossed ("fan") series 
of diagrams. In general they are functions of time and 
momentum. However, in the RMT, or ergodic, limit the 
main contribution to observables is done by the zero mo- 
mentum mode. In this limit the diffuson D^g it, t'; rf) and 



the cooperon C"^(?7, ry'; t) are given by: 



I?;f(i,t';ry) = (1/2) i?(t, i'; ry) x 






c;^iv,v';t) = Cif^,v';t)x 



(7) 



(8) 



{aj\Pom 



2(77-77') 
e ^^^ {a-f\Pi\pS) 



where Po,i are the projectors on the subspaces with the 
total spin S' = 0, 1, respectively: 

{aj\Pom = al^spn, (9) 

{aj\Pim = {Sa,6sp + aZ^a^sp + <^^h) /2. (10) 
D{t,t';T]) and C{ri,i]';t) are defined as follows: 

D{t,t';r]) ^ 9{t-t')exp\- f -f{t" ,Tj)dt"\ , (11) 



C{fj,r^';t) = %-r/)exp|-i^''7(t,r/')dry"|(12) 

where j{t, rj) is determined by the external field: 

7(t, rj) = rm + 7^/2) - 0(t - ,7/2)]2. (13) 

Note that Eqs. (|7I8(I retain crossover triplet terms. 
At times larger than Tso they decay exponentially, and 
the results reduce to those for the RMT model |5J) with 
Hq from GSE, which corresponds to t^o — > 0, so that 
in this model the triplet term would be absent from the 
very beginning. The RMT model (01 with Hq from GOE 
and V from GSE would correspond to a finite Tso- 

The diagrammatic technique in the Keldysh represen- 
tation allows to calculate the time- and energy-dependent 
electron distribution function /(e,t), from which one de- 
duces the energy absorption rate: 



W{t) 



d_ 
dt 



e/(e)de. 



(14) 



The expansion in the number of the diffuson or cooperon 
loops corresponds to the expansion in the powers of the 
mean level spacing S, which is assumed to be the smallest 
energy scale of our problem |8|, [ij • The leading contri- 
bution is given by diagrams containing no diffuson or 
cooperon loops; it corresponds to the Ohmic absorption 
with the rate 



Wo 



2r 



idt<i^y 



(15) 



where the overline denotes the average over the period. 
The next order correction to the Ohmic absorption rate 
is obtained by taking into account one- loop diagrams: 

W{t) = Wo + - f (Se^ " ^) ^ 

X 9t0(i) dt^t ~ 77) C(77, -77; t - r^/2) d77.(16) 

For 4){t) — cosLot the long-time behavior {t ^ 1/r) of 
the above expression takes the following form: 



Wo 



= 1 



{l/2)^/t/U, Tso <i< t. 






2J2- 

(17) 



These two limiting cases differ by the presence or absence 
of the triplet contribution and correspond to the orthog- 
onal or symplectic symmetry classes. Thus the weak lo- 
calization correction to the classical absorption rate in 
the symplectic case is positive and its magnitude is half 
that for the orthogonal case. Exactly the same holds 
for the weak localization correction to the conductivity 
of a quasi- Id disordered wire [iflj, so the relation is 
extended also to the symplectic case. 

The perturbative result (fTTI is valid only for times 
i <C i*. At i » i* the relation ^ suggests to use the 
known results for a quasi-ld wire 111, where, in spite of 
the positive sign of the first weak localization correction, 
all states are localized, with the localization length four 
times larger than in the orthogonal case. Localization 
implies D{u!) -^ at u; ^ 0, for which Eq. ||2J) gives 
W{t) ^ at t ^ cx). 

Quantum kicked rotor with spin. — The exact time de- 
pendence of the correction H17|l is derived rigorously for 
the disordered model Q). However, one can expect that 
the increase of the energy absorption rate in time is a 
general result valid for an ergodic dynamical system pos- 
sessing symplectic symmetry. In order to support this 
statement we introduce a spin degree of freedom into the 
standard kicked rotor model 'l8'|. Correspondingly, V{9) 
in the Hamiltonian Q) acquires a 2 x 2 matrix structure: 

V{6) ~ cos a cos /3 cos 6* -I- 

-I cos a sin/? sin 26* + cr^ sin a sin 0, (18) 

where the parameters a and /3 allow to switch between 
different symmetry classes. Note that the presence of 
both a^ and a^ terms with different dependencies on 6 is 
essential, as otherwise the spin sector of the Hamiltonian 
could be diagonalized by a global spin rotation. 

In order to observe the dynamical anti-localization 
we fix the parameters a = 0.187, /3 — 1.284 which 
corresponds to the symplectic symmetry class [l^l and 
study the evolution of the wave-packets in the momen- 
tum space. The other parameters are T — 1, I — 
10(\/5 - l)/(27r), K = 10/, and the Hilbert space size 
(^max = 16384. This model is very convenient for numeri- 
cal study since the application of the Floquet operator to 
a state can be performed by using the Fast Fourier Trans- 
form algorithm as for the standard QKR model. Fig. Q 
shows that the energy absorption rate indeed initially in- 
creases in time until the strong localization changes this 
behavior to the opposite one. For comparison we plot the 
energy absorption rate for the standard QKR calculated 
for the same values of T, I and K as well. Here, in con- 
trast, the quantum correction is negative from the very 
beginning. 

Conclusions. — The main results of the paper are rep- 
resented by Eq. H17|l and Fig. El The former shows 
that the analogy between the energy absorption by an 
ac driven chaotic quantum dot and the propagation in 
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FIG. 2: The energy absorption rate as a function of time for 
the standard and symplectic kiclced rotor. The energy E{t) is 
averaged over 100 trajectories with different initial conditions. 



a quasi-ld disordered wire is valid also in the presence 
of the spin-orbit interaction, i. e. for systems of the 
symplectic symmetry class, at least at the level of the 
first weak anti-localization correction. Numerical results 
for the quantum kicked rotor exhibit qualitatively the 
same behavior: weak dynamic anti-localization at shorter 
times, and strong dynamic localization at longer times. 
At the weak anti-localization stage the energy for both 
chaotic systems exhibits a peculiar super-Ohmic growth 
in time, in spite of the fact that no resonances are present 
in the considered systems. 

We thank V. I. Fal'ko for fruitful discussions. 
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